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ABSTRACT 


A computer simulation of the ionospheric experiment of the Apollo-Soyuz Test 
Project (ASTP) has been performed. ASTP is the first example of USA/USSR coopera- 
tion in space and is scheduled for summer 1975. The experiment consists of per- 
forming dual-frequency doppler measurements (at 162 and 324 MHz) between the 
Apollo Command Service Module (CSM) and the ASTP Docking Module (DM), both 
orbiting at 221-km height and at a relative distance of 300 km. A network of ground 
stations will also operate and will collect differential and rotating doppler data in DM- 
to-ground radio paths. The computer simulation has shown that, with the doppler 
measurement resolution of approximately 3 mHz provided by the instrumentation (in 
10-sec integration time), ionospheric-induced doppler shifts will be measurable 
accurately at all times, with some rare exceptions occurring when the radio path 
crosses regions of minimum ionospheric density. The computer simulation has 
evaluated the ability of the experiment to measure changes of columnar electron 
content between CSM and DM (from which horizontal gradients of electron density at 
221-km height can be obtained) and to measure variations in DM -to-ground columnar 
content (from which an averaged columnar content and the electron density at the DM 
can be deduced, -under some simplifying assumptions). The simulation has confirmed 
the expectation that simultaneous measurements of horizontal gradients and space-to- 
grounl col umn ar content substantially increase the accuracy of data inversion. 
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DOPPLER MEASUREMENTS OF THE IONOSPHERE ON THE OCCASION 
OF THE APOLLO-SOYUZ TEST PROJECT 

Part I: Computer Simulation of Ionospheric -Induced Doppler Shift 


Mario D. Grossi and Ray H. Gay 


1. INTRODUCTION 

A spacecraft-to-spacecraft doppler-tracking experiment will be performed by the 
Smithsonian Astrophysical Observatory on occasion of the forthcoming Apollo-Soyuz 
Test Project (ASTP). The ASTP mission, scheduled for summer 1975, is the first 
example of USA/USSR cooperation in space. The primary goal of the doppler-tracking 
experiment is to refine the anomalies of the earth’s gravity field. 

The experiment consists of measuring, with a doppler-tracking radio link, the 
relative velocity between the ASTP docking module (DM) and the Apollo command 
service module (CSM). The DM and CSM will both orbit at 221-km height and will 
stay at a relative distance of 300 km for approximately 24 hours (the duration of data 
taking). From doppler signals in a dual -frequency VHF link from the DM to the CSM, 
the relative -velocity data will be inverted into anomalies of the earth’s gravity field 
with a threshold sensitivity of the order of 10 mgal. The integration time for the 
doppler measurements will be 10 sec. 


This report covers the computer simulation of the ionospheric-induced doppler shifts 
(Part I). Part R (SAO Special Report No. 367) covers the inversion of differential 
and rotating doppler shifts. 

This research was supported by NASA under contract NA§_9- 13837 to the Smithsonian 
Astrophysical Observatory. * ^ ^ 
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This experiment, known by the code name ASTP/MA-089, is the first attempt at 
testing the so-called '’low-low" satellite-to-satellite doppler-tracking approach for 
gathering data on medium-wavelength features (horizontal wavelength from about 300 
to approximately 1000 Ion) of the earth’s gravity field. The method, called low-low to 
emphasize the low height of the two orbiting terminals of the doppler-tracking link, 
has been described in the literature by Wolf (1969), Comfort (1973), and Schwartz (1970). 
Its comparative merits with respect to the so-called low-high satellite-to-satellite 
doppler-tracking scheme [described in the literature by Von Bun (1972) and named 
this way because one terminal is in a high orbit and the second in a low orbit] are still 

a subject of debate. 

The Apollo-Soyuz Test Project will give us the first opportunity to compare the 
two approaches. In fact, in addition to the low-low method tested with the MA-089 
experiment, Goddard Space Flight Center will be conducting an on-board ASTP low- 
high experiment, known under the code name of MA-128. 

hi the MA-089 experiment, the doppler error induced in the measurements by the 
ionosphere will be corrected by a dual-frequency method (Guier and Weiffenbach, 1960). 
The spacecraft-to- spacecraft doppler-tracking link will, in fact, operate at two coherent, 
harmonically related frequencies, 162 and 324 MHz, and the effect of the ionosphere 
will be quantitatively determined by exploiting the phenomenon of ionospheric frequency 
dispersion. These frequency values have been chosen because instrumentation has 
already been developed and is available and ground stations tuned at 162 and 324 MHz 
have already been deployed. 

By means of a transmitter installed on board the DM as a dual -frequency beacon, 
an ionospheric experiment (Grossi, 1974) will also be performed (see Figure 1), in 
which differential doppler data wall be collected in DM-to-CSM paths and differential 
and rotating doppler data will be collected in DM-to-earth paths. The doppler meas- 
urement resolution will be of the order of 3 mHz, with a 10- to 100 -sec integration 

time. 

The use of the DM as a terminal of the MA-089 doppler-tracking link was initially 
suggested by NASA Astronaut Commander Eugene Cernan. 
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Although these ionospheric measurements are secondary to the gravity-field 
measurements of the MA-089 experiment, chronologically they were proposed to the 
National Aeronautics and Space Administration earlier (Grossi, 1972), when the ASTP 
mission was known as the International Rendezvous and Docking Mission 1975 
(IRDM 1975). 
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2. SCIENTIFIC GOALS 


The MA-089 ionospheric experiment represents an extension of well-known methods 
(differential doppler and rotating doppler) of measuring the integrated electron concen- 
tration and other properties of the ionosphere along a radio path between a terminal 
moving in the ionosphere or above it and a station on the ground. Even before 
artificial-satellite flights, these techniques wore in use in suborbital rocket flights 
(Seddon, 1953; Jackson and Seddon, 1958) performed to measure ionospheric param- 
eters. Since 1957, a wealth of literature has appeared, based on the use of multifre- 
quency doppler links between satellites (orbiting inside the ionosphere and above it) 
and ground stations (Al’pert, 1958, 1965; Garriott, 1960a, b; de Mendoza, 1963; Misyura 
et al. , 1964). Recent additions to the literature include papers byAl'pert (1973, 1974) 
and Tyagy (1974). 

The fundamental problem that users of space-to-ground doppler links must face 
is the fact that the differential doppler shift observed when receiving two coherent, 
harmonically related frequencies is directly connected to the time derivative of the 
columnar electron content and not to the columnar content itself. This difficulty can 
be removed if simultaneous measurements of the Faraday rotation (also called rotating 
doppler shift) affecting the link are carried out. men the rotating doppler shift is not 
observed, the problem is under deter mined, and the accuracy obtainable from inverting 
the differential doppler data into columnar content strongly depends on the presence of 
horizontal gradients in that region of the ionosphere that is swept through by the space- 
to-ground radio path while the space-borne terminal is in motion. 

The MA-089 experiment introduces a new measurement feature, which is connected 
to the problem of horizontal gradients: namely, the measurement of their values at a 
height of 221 km (the common orbital height of both the DM and the CSM) in the orbital 
plane of the DM/CSM pair. These measurements are performed simultaneously with 
space-to-ground differential doppler measurements, by using the DM as a dual- 
frequency beacon. 



As will be shown in Part II of this report, a significant increase in accuracy in 
the inversion of differential doppler data into columnar content is achievable. 

In summary, the MA-089 experiment will yield the following: 

A. Changes of the columnar electron content between the two spacecraft, from 
which the horizontal gradient of the electron density at a height of 221 km (along the 
orbital path of the DM/CSM pair) will be obtained. 

B. Variations in the ground-to-spacecraft columnar content, from which an 
averaged columnar content and the electron density at the DM location can be deduced 
under some simplifying assumptions. 

In addition, it will be possible to try to probe traveling ionospheric disturbances 
(TID) with both the space-to-ground link and the DM-to-CSM link. This will enable us 
to measure the TID electron-density fluctuations with respect to the surrounding 
background, to sample the TID wave structure, and to monitor the motion of the TID. 

Finally, the experiment will attempt to detect the boundaries of turbulent regions 
of the ionosphere, such as the aurora oval and the equatorial sporadic E belt. 
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3. MEASUREMENT APPROACH 


3. 1 Description of the Method 

Let us consider first a simple case: the vertical ascent in the ionosphere of a 
spacecraft that carries the transmitter (with pulsation of a doppler-tracking link. 
The receiving terminal is on the ground. When the spacecraft is at a height r(t), the 
link’s phase-path length is 

r(t) 
co, r 

*l (t) = c J n l (r) dr ’ 

0 

and its time derivative is 

CO r W 

*l«-T 3t / n l <r)dr ’ (2) 

0 

where c is the velocity of light in free space and n(r) is the index of refraction at 
height r. The change of phase path with time is due in part to the spacecraft motion 
and in part to temporal changes of the index of refraction along the vertical between 
the spacecraft and the ground. If there are no temporal changes along the path, then 
by assuming the spacecraft to be at height r Q at time t Q , we have 

*i (t o ) = ~c • (3) 


If there are temporal changes, equation (2) becomes 


co i 


*1<V = if 


r(t) 


9n, (r) 


.at — dr 


W 
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r (t) 

where the term J [8n^(r)/9t] dr represents the temporal variation of the columnar 
0 

refractivity in the entire vertical path between the spacecraft and the ground. The 
problem is underdetermined, and unless this variation is otherwise measured or 
becomes negligibly small, the inferring of n^r^) from $ ^ (t Q ) is affected by error. 

When the effects of the earth's magnetic field and the collision frequency on the 
index of refraction are disregarded, equation (3) can be rewritten as follows: 


CO. 


= C r <V 



(5) 


where N is the electron density and e and m are the charge and the mass of the electron, 
respectively; 2-n-e /m = 1587. 6, when N is in electrons m and is in rad sec . 

From equation (3), we can determine the local index of refraction (and, hence, the 
electron density) at the spacecraft height by monitoring the received doppler shift and 
by knowing, independently, the ascension velocity of the spacecraft and the frequency 
radiated. When operating with a single frequency, it is necessary to know these two 
parameters very accurately. Since the ionospheric -gene rated doppler shift is a few 
millihertz, we must know the spacecraft velocity with an accuracy better than a fraction 
of a millimeter per second (when the link operates with a wavelength of, say, 1 m) and 
the frequency of the link within a fraction of a millihertz (with integration times of the 
order of 10 to 100 sec, the duration of the ionospheric crossing). However, by adding 
a second frequency, oo 2 , that is phase coherent and harmonically related to w , the 
equation of the differential doppler shift in the spacecraft-to-ground link, when tem- 
poral changes of the columnar refractivity between the spacecraft and the ground are 
neglected, becomes 


00 n 


r (t) 


6$ = $ x (t) - 


~ *,<t) = 

co 2 2 


= -r at J [n i< r > - n 2 < r >i dr = -r *<V [n i< r o> - W 1 


( 6 ) 
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If we disregard the refractive effects of the earth's magnetic field and the collision 
frequency, we get 




2ue 2 N 

m 


2 2 

“l " f °2 

2 2 

0), to., 


CO 


From equation (7), we can now deduce that the contribution to the measurement error, 
arising from an error in estimating the link's frequencies, is virtually eliminated. 

We see from equation (5) that an error in frequency contributes directly to the 
error of (t^). However, from equation (7), we find that the error in frequency 
must now be multiplied by the quantity to 2 - 1) Uj 2 , which is usually very small. For 
instance, if « /« 2 = n = 2 and ^ X 300 X 10 8 rad sec" 1 , we have 
(n 2 - 1) oj^ 2 = 8. 6 X 10" 19 , and therefore the influence of the frequency error in the 

overall error of 4^) is eliminated. 


When the spacecraft trajectory is no longer vertical, other analytical expressions 
will have to be employed. Alternative formulas have already been developed by 
several authors and are available from the literature. For the case in which horizontal 
gradients of the electron content are negligible, we can write (following Al'pert, 1973): 


64 (t) - a Q 


cos <j> s 
Y s 


N + r +• 


s COS 4> 


s 1 N 


( 8 ) 


where the velocity components correspond to the height z g of the transmitter, N g is the 
local value of the electron concentration at that height, 


Mf 


N dz 


(9) 


2 

is the mean value of the integrated electron concentration in a column of 1-cm cross 
section, and the coefficient a Q is 
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In equation (8), r(t) denotes the time-dependent radius vector joining the point of 
observation to the moving source, which is assumed to be approaching the observer 
and located at height z(t). Moreover, the radial, horizontal, and vertical velocity 
components of the source are denoted, respectively, by r(t) (along the line of sight), 
x(t), y(t), and z(t). 

In the general case, we can write for the electron concentration 

N = N(R,0,x,t) = N(z,x,y,t) , < n > 


where (see Figure 2) 

R = R Q + z , x = R q 0 , y = R q X > < 12 > 

and (z, 0, x) is a variable coordinate in a spherical orthogonal coordinate system along 

the wave -propagation trajectory joining the point of observation (0, 0, 0) and the point 

fz 0 ,\) where the transmitter is located at time t. The difference in the doppler 
' s’ s’ s 

frequency shifts for the two coherent radio waves can then be written in the form 


R s (t) 


siw =a o 5 J • 


(13) 


R, 


If certain conditions are satisfied, as happens in the case of the MA-089 experi- 
ment, we obtain (Al’pert, 1973) 


- N S S5TT + N x -Vs' N t 


6$ = a, 


(14) 




Figure 2. Nomenclature used in equations (11) through (15). 
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When the further simplifying assumption is made that the medium is plane parallel 
(i.e., the earth’s sphericity is neglected), the various parameters of equation (14) 
are defined as follows (Al’pert, 1973): 


N = — f N dz + t- f P- z dz , 

x z J z sin 4> n cos <{> n J 9x 1 

S ,rv S U " a 


N 

y z g cos Y() 


g 

i r 

os (j> n J 


9N 

ay 


z dz , 


N = -At- ffdz . (15) 

t cos 9 ^ J 9t 

Even in the planar approximation, the problem is underdetermined. This difficulty 
is, in part, alleviated in the MA-089 experiment by the fact that the DM-to-CSM dual- 
frequency link measures the quantity 9N/8x (the horizontal gradient in the orbital plane) 
at the ASTP orbital height: 


CO , ■ CO 

6 *=*i<t)--*2 = -; 


CSM 


CSM 


9_ 

at | 


J" n 1 (x) dx - J* n, 

DM DM 


(x) dx 


(16) 


Assuming the two spacecraft remain at a constant separation (relative velocity = 0), 
we obtain 


CSM 

6 * = -£ / b [n i (x) - n 2 (x>1 If ** • (17) 

DM 

By neglecting the refractive effects of the. earth’s magnetic field and the collision fre- 
quency, we can rewrite equation (17) as follows: 
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(18) 


8$ = — 


o/2 2 

“l 2ne 2 /“I - “2 


CSM 


m 


2 2 
W 1 W 2 


r m 

J ax 


9N 9x 

9x at 


dx 


DM 


where 9x/8t is known from orbital -mechanics considerations and 9N/8x can therefore 

be obtained from the measured values of 6$ in the DM-to-CSM path. Before it can be 

used in equation (15), however, 9N/9x must be known all along the vertical z. hi 

/*z 

fact, what is needed is the function I (9N/ax)(z) dz (and not just 9N/9x at the ASTP 
orbital height of 221km). We must therefore construct a model of 9N/9x (z) in the 
lower ionosphere, with the constraint of satisfying both the value measured at 221 km 
by the DM-to-CSM link and a value equal to zero at the ionosphere's bottom. A linear 
variation of the gradient between these two values thus seems to be an acceptable 
assumption. 


In the DM-to-CSM link, the transmitting antenna is linearly polarized and the 
receiving antenna is circularly polarized. Therefore, we cannot observe the Faraday 
rotation (rotating doppler) phenomenon. In the DM-to-ground link, though, we plan 
to adopt on the ground a receiving antenna that is linearly polarized and expect to be 
able to monitor the amplitude of the received signal. If these expectations materialize, 
we can observe the quantity 6$ H , which, for quasi-longitudinal propagation (Tzedilina, 
1962; Al’pert, 1965), can be expressed 




moo 



(cos 9) N 


dz 

COS cj> 


9 


(19) 


where 6<j> T . is the frequency at which the amplitude of the received signal changes. In 
H 

the absence of horizontal gradients in the electron-density profile, equation (19) reduces 
to 



2Tre 3 H 

V- tj 

2 2 2 1 
m w c 


cos 0„ z 
s s 

cos <j> 

T s 


N + sin y 
s ‘s 



+ 


cos <j> 


T sin y s 


N 


) 

I 


( 20 ) 
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where y is the angle between H Q and the vertical z, £ is the angle between the projection 
of Hq along the y axis in the yz plane and the xz plane, and 

cos 0 = cos <{> cos y + sin <j> sin y cos £ 

We can now obtain N and N at any given point by means of equations (14) and (20). 

s 

If we no longer assume the absence of horizontal gradients in the electron-density 
profile, we must go back to the general formula (19). 

3. 2 Evaluation of the Measur e ment Error 

3.2. 1 Instrumentation contribution 

The instrumentation error in the MA-089 experiment is dominated by the following 
two sources of inaccuracy: 

A. Oscillator instabilities. The frequency- stability curve in Figure 3 indicates 

12 

that each oscillator has a stability of approximately 1. 5 parts in 10 (for 10- to 100- 
sec integration time). Because each link of the doppler-tracking system uses two 
oscillators, we have, at 324 MHz, an error 

Af = 324 X 10 6 X 1. 5 X 10~ 12 V2 = 0. 687 mHz. 

At 162 MHz, Af = 0. 3435 mHz. 

B. Phase -locked-loop tracking error. The signal-to-noise ratios (SNR) expected 
to be available in the doppler-tracking link are 

~32 db at 162 MHz , 

~26 db at 324 MHz . 

The. tracking-error variance in the phase-locked loop, therefore, will be 

a 2 = 6. 3 X 10“ 4 rad 2 at 162 MHz , 

P . 

a 2 = 2. 5 X 10" 3 rad 2 at 324 MHz . 

P 
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FREQUENCY STABILITY 



Figure 3. Curve of the stability of the oscillators used in the ASTP doppler-tracking 
experiment. 
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Since the lower frequency is doubled before it is sent to the doppler processor, the 
phase -error variances for the two channels are identical; and since the loop bandwidths 
are both 5 Hz, the phase errors observed at time intervals T. = 10 sec are independent. 
The total root-mean-square error, therefore, is 


V2 

CT f °p 2ttT. 


V2X, 10~ 2 X 2.5 
2ir X 10 


= 0. 565 mHz 


The total instrument error in taking raw doppler data at 324 MHz is thus approximately 
0. 8 mHz. 

We performed laboratory tests on the ASTP doppler-tracking instrumentation with 
the transmitting and receiving terminals connected back to back by means of an atten- 
uator that simulates path losses. The following errors were typically measured; 

2. 04 mHz at 162 MHz , 

1. 89 mHz at 324 MHz , 

1. 25 mHz for the differential doppler shift . 

The increase in error above the expected value is attributed to the inability of the 
oscillators to perform as well as they could in principle. This is not, however, the 
sole reason for downgrading our measuring-accuracy ambitions. As we show in Section 
3. 2. 2, even with oscillators performing as expected, multipath propagation effects 
would cause the doppler error vastly to exceed the value of 0. 8 mHz, initially estimated 
to be the overall error for the instrumentation. 


3.2.2 Multipath propagation effects 

Because the antennas used at the two terminals of the doppler-tracking link have 
low directivity (the transmitting antenna on the DM has an average gain of -5 db, and 
the receiving antenna on the CSM, a gain of +6 db), reflections from the earth's surface, 
as well as the DM-to-CSM direct ray, will reach the CSM antenna. This is a well- 
known multipath phenomenon. 
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A simplified analysis of this effect can be carried out by assuming that the multi- 
path structure contains only two rays (Figure 4) and that reflection on the earth’s 
surface takes place at the midpoint R. 


Because the electric vector of the transmitted signal is perpendicular to the 
orbital plane of the spacecraft, the direct wave and the wave incident on the earth's 
surface are horizontally polarized. If the electric field at the CSM (see Figure 4) 
arriving from the DM along the direct path l ^ is E^, the reflected electric field at 
CSM arriving along the indirect path is 


E 


r 



DH h E 0 e 




) 


where D, the defocusing factor due to the earth's curvature, is given by 


D = 


(* + 2r sin e/j) (* + 2r ). 


- 1/2 


and ^ are the magnitude and phase of the Fresnel reflection coefficient for hori- 
zontal polarization; and 0 is the phase-path difference, given by 


9 =¥< S -’l2> ' 


The total field at the CSM is 


12 


j<V 0 > 


E, , = E + DR u E»e , 

tot 0 s h 0 ’ 


from which it follows that 


E tot! 2 ( l 12 \ 2 *12 

-“t- 1+ -r DR h + 2 -r D \ cos (f-h - 9) , 


E 


O 1 


17 



DM 


CSM 



Figure 4. Multipath configuration. 
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where 


1,2 + J 

r2 + r 


2i n r 

2 




1/2 1 1/2 


2 , 2 

Assuming isotropic antennas at both ends of the link, we have plotted ( /|Eq| 
in Figures 5 through 8 for desert and sea water, for two frequencies. Values of 
| ^|/| Eq j are also given in column 4 of Tables 1 thi'ough 4. It follows from the 

above equations that the normalized reflected field at the CSM is 




Col umn 5 of Tables 1 through 4 lists some values of |E^|/|E 0 | . The tables and 
graphs indicate that the worst situations occur when reflection takes place on sea 
water. It can be shown that the worst-case SNR degradations (see Figures 7 and 8 for 
the case of 300-km range and 56° grazing angle) are 

-8 db at 324 MHz , 

-6 db at 162 MHz . 


The multipath problem can be alleviated by tilting the CSM in such a way that, in 
the direction from which the ray reflected from the ground is received, the antenna 
has a substantially smaller gain than in the DM/CSM line of sight. Figures 9 through 
12 indicate that when the CSM antenna has a 10-db gain difference in the two directions, 
the degradation due to sea-water reflection is less than 1 db at 162 and 324 MHz. For 
a 15-db gain difference, the multipath effect on the link' s SNR becomes totally 
negligible. 

The problems created by the multipath phenomenon are not, however, limited 
solely to the deterioration of the SNR, which is largely limited and therefore inconse- 
quential. Unfortunately, the multipath problem also unfavorably affects the behavior of 
the phase-locked loops. According to an analysis performed by Stiff ler et al. (.1975), 
the error due to this last effect can be evaluated as 2.7 mHz (in 10-sec integration time) 
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RELRTIVE DIELECTRIC CONSTRNTs 3.00 EO 
SURFACE CONOUCTIVITY= 1.00 E -6 CRHOS/M) 
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Figure 5. Signal fluctuations at the receiver terminal (162 MHz, DM-to-CSM link above desert). 
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Figure 6. Signal fluctuations at the receiver terminal (324 MHz, DM-tc-CSM link above desert). 
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Figure 8. Signal fluctuations at the receiver terminal (324 MHz, DM-to-CSM link above sea water) 
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Table 1. Multipath characterization (162-MHz link above desert) 


to 


Surface dielectric constant = 3 

-5 

Surface conductivity =10 mhos m 

Frequency = 162 MHz 


Height of both spacecraft =221 km 

Radius of earth = 6370 km 

Maximum separation of spacecraft = 500 km 


Grazing angle 

Direct 
path length 
(km) 

Indirect 
path length 
(km) 

l E totl 

1 E r 1 

i E ol 

(db) 

K\ 

(db) 

83° 50 

50.35 

444. 76 

-0.02 

-30. 60 

77.00 

101.95 

453.22 

0.07 

-24. 45 

71. 00 

151.89 

466.54 

-0.80 

-20. 96 

65. 50 

200.73 

484. 06 

-0.72 

-18.51 

60.50 

248.74 

505. 14 

1.19 

-16.63 

55. 50 

301.41 

532. 14 

1.42 

-14. 95 

51.50 

347.93 

558. 91 

-1.87 

-13. 72 

47.50 

399. 47 

591.29 

0.61 

-12. 56 

44.00 

449. 76 

625.24 

-2.46 

-11.59 

41.00 

497.65 

659.40 

-1.01 

-10.79 

- 

499. 99 

661.11 

- 
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Table 2. Multipath characterization (324-MHz link above desert) 


to 

CJI 


Surface dielectric constant = 3 

-5 

Surface conductivity = 10 mhos m 

Frequency = 324 MHz 


Height of both spacecraft =221 km 

Radius of earth = 6370 km 

Maximum separation of spacecraft = 500 km 


Grazing angle 

Direct 
path length 
(km) 

Indirect 
path length 
(km) 

l E totl 

|e j 

i r i 

"TV 

<db) 

TV 

(db> 

83°. 50 

50.35 

444. 76 

0.25 

-30. 60 

77.00 

101.95 

453. 22 

0.50 

-24. 45 

71.00 

151.89 

‘466.54 

-0.75 

-20. 96 

65. 50 

200.73 

484. 06 

0.08 

-18.51 

60.50 

248.74 

505. 14 

-1.35 

-16.63 

55.50 

301.41 

532. 14 

-1.65 

-14. 95 

51.50 

347. 93 

558.91 

-1.49 

-13.72 

47.50 

399. 47 

591.29 

1.73 

-12.56 

44.00 

449.76 

625.24 

-1.94 

-11.59 

41.00 

497.65 

659.40 

1.25 

-10.79 

- 

499. 99 

661.11 





Table 3. Multipath characterization (162-MHz link above sea water) 


t>0 

05 


Surface dielectric constant =80 
Surface conductivity = 4 mhos m 

Frequency = 162 MHz 


Height of both spacecraft =221 km 

Radius of earth = 6370 km 

Maximum separation of spacecraft = 500 km 


Grazing angle 

Direct 
path length 
(km) 

Indirect 
path length 
(km) 

l E tJ 

|E 0 I 

<db) 

|E r l 

|E 0 I 

(db) 

83? 50 

50.35 

444. 76 

-0.08 

-19. 84 

77.00 

101.95 

453. 22 

0.43 

-13.88 

71.00 

151.89 

466.54 

-2.86 

-10.66 

65.50 

200.73 

484. 06 

-2.32 

-8.56 

60.50 

248.74 

505. 14 

3. 17 

-7.06 

55.50 

301.41 

532. 14 

3.55 

-5.85 

51.50 

347. 93 

558.91 

-5.74 

-5. 03 

47.50 

399. 47 

591.29 

1.93 

-4. 33 

44. 00 

449. 76 

625. 24 

-7.67 

-3. 80 

41.00 

497.65 

659. 40 

-1. 60 

-3.41 

— 

499.99 

661.11 

— 

— 



Table 4. Multipath characterization (324-MHz link above sea water). 


Surface dielectric constant = 
Surface conductivity = 

Frequency = 

80 

4 mhos m * 
324 MHz 

Height of both spacecraft 
Radius of earth 

Maximum separation of" spacecraft 

Grazing angle 

Direct 
path length 
(km) 

Indirect 
path length 
(km) 

l E totl 

Kl 

<db). 

|E r l 

l E ol 

(db) 

83°. 50 

50. 35 

444. 76 

0.80 

-20. 14 

77.00 

101.95 

453.22 

1.51 

-14. 17 

71.00 

151.89 

•466.54 

-2.44 

-10.94 

65. 50 

200.73 

484.06 

2.56 

-9. 17 

60.50 

248. 74 

505. 14 

-4. 48 

-7. 32 

55.50 

301.41 

532. 14 

-5.27 

-6.09 

51.50 

347. 93 

558.91 

-3.40 

-5. 26 

47.50 

399. 47 

591.29 

3.92 

-4. 55 

44.00 

449. 76 

625.24 

-4.80 

-4.01 

41.00 

497.65 

659. 40 

2.90 

-3.60 

_ 

499. 99 

661. 11 

- 
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Figure 9. Signal fluctuations at the receiver terminal (162 MHz, DM-to-CSM link above desert, 10-db attenua- 
tion in the direction of the reflected path). 
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Figure 10. Signal fluctuations at the receiver terminal (324 MHz, 
uation in the direction of the reflected path). 
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Figure 11. Signal fluctuations at the receiver terminal (162 MHz, DM-to-CSM link above sea water, 10-db atten 
uation in the direction of the reflected path). 
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Figure 12. Signal fluctuations at the receiver terminal (324 MHz, DM-to-CSM link above sea water, 
attenuation in the direction of the reflected path). 
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when the difference in antenna gain between the direction of the direct path and that of 
the reflected path is 10 db, and 1. 35 mHz when the gain difference is 15 db. 

3. 2. 3 Overall expected error 

From the above discussion, we conclude that the overall expected error in meas- 
uring both the doppler shift at a single frequency and the differential doppler shift is 
between 2 and 3 mHz. We shall adopt 3 mHz as the error in the raw data. 

The estimate of the overall measurement error in the complete MA-089 experi- 
ment must include the contribution to the error budget of data processing and data 
inversion. This estimate will be performed in Part II of this report. 
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4. DOPPLER-TRACKING INSTRUMENTATION 


Figures 13, 14, and 15 show block diagrams (Stiffler et al . , 1975) of the beacon 
transmitter installed on the DM (0. 1-watt output power at each of the two unmodulated 
carriers) and of the two-frequency receiver, doppler processor, and recorder, all 
installed on board the CSM. The antenna on the DM is a linear monopole oriented 
along the axis of maximum moment of inertia of the DM, perpendicular to the orbital 
plane. As initially suggested by Colombo (1974), a maneuver of the CSM/DM (before 
separation) will impose on the DM a rotation of approximately 3° sec 1 around that 
axis. It is expected that this orientation will be maintained by the monopole antenna 
during the lifetime of the experiment. The receiving antenna on the CSM is a circularly 
polarized microstrip panel structure. 

On board the CSM, the 10-sec doppler samples are encoded in digital words 
organized in 73-word frames, each of which contains a frame-number word. Each 
frame contains a total of 5063 bits and represents approximately 12 min of data. There 
will be, therefore, 120 frames in 24 hours of data collection. The frames will be 
stored in a redundant pair of on-board, single-channel, start-stop tape recorders, 
where the elementary bits are written in with a speed of approximately 1 kbit sec 1 . 
Every frame will utilize approximately 5 sec of tape plus a guard band of some 7 sec. 
Thus, a full 24 hours of data collection will be compressed into approximately 25 min 
of actual tape recordings. Some similar arrangement is foreseen for the ground 
recordings. 

Because the doppler processor is the heart of the doppler-tracking system, it will 
be described hereunder in more detail. The concept of the processor has been worked 
out by Weiffenbach (1973). 

The processor consists of five counters: a 10-sec interval counter, two up-down 
counters, and two zero-crossing counters, with one up-down counter and one zero- 
crossing counter associated with each of the two channels. The operation of the proc- 
essor can be explained by referring to Figure 16. The 10-sec counter identifies the 
time instants t Q , t 0 = t Q + 9. 996 sec, t^ = t Q + 10 sec, 10 sec, etc. The points 
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Figure 14, Block diagram of the receiver, 
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Figure 15. Block diagram of the doppler processor 
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of and t g are determined by the first positive-going zero crossing in the heterodyned 
signal occurring after times t^ and tg, respectively. (Note that while t^ and tg are 
common to both channels, there are two t^s and two tg’s, one of each associated with 
each of the two channels.) For each channel, the associated up-down counter is 
enabled at time t Q and counts the number of microsecond clock pulses observed before 
it is disabled at time t^ When it is enabled again at t g , this time in the down-count 
mode, it down-counts the number of clock pulses in the interval (t^, t^). Simultaneously, 
the zero-crossing counters count the number of positive -going zero crossings in the 
interval (t , t ) for each channel. These time and zero-crossing counts constitute the 
raw data associated with the observation interval in question. The counters are reset 
and the whole process begins anew at time t^. 

Four quantities, namely, r F = [t ^ (i) - t^] - [tg(i) - tg] and Z^, the number of 

positive -going zero crossings in the interval (t^, tg), are thus generated during 

each 10-sec observation period. (Here, i = 1, 2 is an index used to distinguish 

between the two channels. ) The average frequencies observed during this interval are 

estimated by f. = Z./(9. 996 - T.); we discuss the error in this estimate below. Although 
ill " 

it would be possible to design the processor itself to determine f., it was found to be 
more efficient to record the data as follows: The first word in each frame contains the 
quantities Z 9 , Z , Tg, and T^ all remaining words contain AZ 2 , AZ^ T 2 , and T^, 
with AZ defined as the difference between the current value of Z and its immediate 
predecessor. 

Thirteen bits each (12 bits plus sign) are used to record T l and Tg, 15 bits each 
for Z and Zg, and six bits each for AZ X and AZg. This allows unambiguous meas- 
urements of all frequencies in the range 250 to 1750 Hz (+7 50 -Hz deviation from the 
1000-Hz center frequency), with frequency variations as large as ±3. 2 Hz between 
observation intervals and up to ±1. 6 Hz from channel to channel. These ranges are 
all well in excess of the maximum anticipated for these quantities. 

Because the integrity of these data is vital to the success of the MA-089 experi- 
ment, various encoding techniques for protecting against random substitution and 
insertion and deletion errors were examined. The use of two tape recorders, of 
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course, provides considerable error protection, but not against errors induced before 
recording, nor does it offer any protection should one of the two recorders fail. We 
decided that a modified, shortened Hamming code would provide significant error pro- 
tection at the cost of a negligible increase in the complexity of the processor and a 
modest increase in its memory requirements. Adjoining seven parity-check bits to 
each word as it is read into the processor’s buffer memory allows any single substitu- 
tion errors to be corrected, any double (and most multiple) substitution errors to be 
detected, and insertions or deletions of up to three bits in any word to be found and 
confined to those words directly affected. 

The first word in each frame (after the eight-bit frame number) is then 87 bits 
long (80 information bits plus seven parity-check bits), with the remaining words 
requiring 69 bits. As already mentioned, a frame contains 73 words, resulting in a 
grand total of 5063 bits. Figures 17 and 18 illustrate the frame and word formats. 

In addition to the doppler information, each data word contains rate-gyro infor- 
mation obtained from the roll, pitch, and yaw spacecraft channels (each with eight-bit 
accuracy). This makes it possible to remove from the data the doppler shift due to 
CSM rotational motions, or at least to identify the times when its attitude jets are 
fired. 

In order to mark in the record the data invalidated by a loss-of-lock condition in 
either of the two tracking loops, the first affected word is recorded with an 01111 pre- 
fix and is followed by the appropriate number of "0's" rather than "l's. ” The next 
word is shortened, and any subsequent invalidated data are recorded in the same 
manner as the first until lock is regained. Data from the first valid observation fol- 
lowing this event are then recorded in the 46-bit format, and the normal format is 
again resumed for the remainder of the frame. 

It is tentatively planned to turn on the DM-borne dual-frequency transmitter on 
July 23, 1975, at the time of the DM/CSM separation. Its battery-limited emissions 
are expected to last approximately 24 hours. Ground receptions will be attempted at a 
few sites, one of which is possibly contained in the orbital plane, for at least one pass 
in the 24 hours. 
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Figure 17. Data-frame format. 
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Figure 18. Data- word format. 
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5. COMPUTER SIMULATION OF IONOSPHERIC-INDUCED DOPPLER SHIFTS 
5. 1 Background 

In order to check the functioning of the software being developed for the reduction 
and processing of the MA-089 experiment data, we have performed a computer simula- 
tion of the expected doppler measurements, anticipating as closely as possible the 
results of the experiment. 

This simulation has been accomplished by using both a model of the 1975 summer 
ionosphere (Odom, 1974) and a ray-tracing computer program (Sforza, 1970). The 
ray-tracing program neglects the effect of the earth's magnetic field but otherwise 
traces rays in a model nonhomogeneous ionosphere, whose structural detail is reason- 
ably realistic. With these tools, we have computed the differential doppler shift 
expected between the DM and the CSM (relative separation = 300 km) while the pair 
moves in a complete orbit around the earth. The orbital plane, inclined -52° with 
respect to the equator, has been chosen to cross the equator at noon and midnight 
local time (LT) and to reach the highest northern latitude at 0600 LT and the highest 
southern latitude at 1800 LT. hi the direct motion of the pair, the CSM has been 
assumed always to be East of the DM, and two ground stations have been assumed to 
be contained in this orbital plane and located at 0430 and 0930 LT. These times corre- 
spond to the lowest and highest ionospheric density as weH as to extremes in horizontal 
gradients. The simulated data we generated will be inverted into the ionospheric 
parameters of the model (see Part II). How close the original and the reconstructed 
data come to each other will indicate the accuracies of the ray-tracing simulation and 
the inversion. Once the actual experiment data become available in late summer 1975, 
we will invert them by means of the same software used to invert the simulated data. 
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5. 2 Model of the Summer 1975 Ionosphere 


Appendix A (Odom, 1974) gives a detailed description of the ionospheric model we 
have adopted. The following quantities appear in the tabulation: 

H t = height (km), 

NL = North Latitude (deg), 

WL = West Longitude (deg), 

_3 

N = electron density (el cm ) 

-3 -1 

DNDR = vertical gradient of the electron density (el cm km ), 

-3 - 1 

DNDNL = latitudinal gradient of the electron density (el cm rad ), 

-3 -1 

DNDWL = longitudinal gradient of the electron density (el cm rad ). 

The model is fully three-dimensional and includes day-night and night-day transitions. 
The highest electron content is along the vertical of ~1000 LT, with a secondary rela- 
tive maximum at 1600 LT. Electron- density minima are along the verticals of 
-0300 and 2300 LT. 

5. 3 Differential Doppler Shift in Spacecraft-to-Spacecraft Paths 

From the ionospheric model described above and in Appendix A, we have used ray 
tracing to compute the differential phase -path length expected at 162 and 324 MHz 
between the DM and the CSM, which are assumed to orbit at the common height of 
221 km and to be separated by 300 km. We define the differential phase-path length 
as the function (in meters) 

S 

$ = J* (|A - 1) ds , 

0 

or the difference between the electrical length of the DM/ CSM separation and its 
geometric equivalent. Because p < 1 in the ionosphere at the ASTP frequency, the 
electrical length (= f p ds) is smaller than the geometric length and $ is therefore 
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negative. Figure 19 provides a plot of this function, with the longitude of the DM as 
the independent variable. Longitude 0 hours is taken to be at local midnight, and the 
adopted position of the orbital plane was given in Section 5.1. Figure 20 shows the 
computed values of the differential doppler shift in the DM-to-CSM path: 

6 $ = $ - 2$ 2 , 

where 4. and <£ 9 are the doppler shifts measured at 324 and 162 MHz, respectively. 

The quantity 5$ can be inverted into horizontal gradients of electron content (see 
Part II). 

5. 4 Differential Doppler Shift in Spacecraft-to-Ground Paths 

Similar calculations have been done for DM -to-ground paths. To cover a variety 
of ionospheric conditions, we have assumed the ground stations to be located in the 
orbital plane and situated at various local times. Figure 21 gives the differential 
phase-path length at 162 and 324 MHz for ground stations at 0430, 0930, and 1230 LT. 
Because the first and third cases represent widely different ionospheric conditions, we 
have computed the differential doppler shifts for each and present the results in 
Figures 22 and 23. 

The space-to-ground differential doppler data will be inverted into columnar 
electron content and electron density at the DM height, as will be done with the actual 
data collected in the summer of 1975 by the MA-089 experiment (see Part II, Sections 3 
and 4). 

5. 5 Modeling Approach for the Rotating Doppler Shift in Spacecraft-to-Ground Paths 

As already mentioned, we expect at least some of the ground stations available 
for the MA-089 experiment to be equipped with linearly polarized antennas. It will 
therefore be possible, by recording the signal's amplitude, to measure the Faraday 
rotation in DM-to-ground paths. 
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DIFFERENTIAL LENGTH <f> ( m ) 



HOURS LONGITUDE 

Figure 19. Differential DM-to-CSM phase-path length. 
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Figure 20. Expected differential doppler shift in the DM-to-CSM path. 











Figure 22. 


EXPECTED DIFFERENTIAL DOPPLER 
IN (DM) * TO-(EARTH) PATHS 



Expected differential doppler shift in DM-to-earth paths , with the ground 
station in the orbital plane at 0430 LT. 
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EXPECTED DIFFERENTIAL DOPPLER IN 
(DM) - TO-(EARTH) PATHS 



Figure 23. Expected differential doppler shift in DM-to-earth paths, with the ground 
station in the orbital plane at 1230 LT. 
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We plan to perform a computer simulation of this phenomenon in the future by 
choosing one of the available ray-tracing programs (Grossi and Langworthy, 1966; 
Jones, 1968) that includes the effect of the earth's magnetic field on the index of 
refraction. The simulation will trace ordinary and extraordinary rays, together with 
their pertinent properties (e.g. , phase-path length and geometric distances). By 
means of this program, we intend to generate simulated values of 6$ H according to 
expressions (19) and (20) of Section 3. 1. 

5. 6 Modeling - Approach for Traveling Ionospheric Disturbances and Ionospheric 
Turbulence 


Methods for detecting ionospheric disturbances from reception of satellite 
emissions have recently been proposed (Bratteng and Forsyth, 1974; Forsyth, 1974). 
The MA-089 experiment offers the opportunity of adding the performance of spacecraft- 
to-spacecraft sounding of the disturbance along roughly horizontal paths. 

A traveling ionospheric disturbance (Thome, 1968; Francis, 1974) could be 
sampled repeatedly at successive passes of the DM/CSM pair. The electron-density 
distribution inside the TID could be measured and some rough tracking of the TID 
position performed. In order to simulate the effect of a TID on the DM-to-CSM doppler 
tracking link, we are planning to add, to our present ionospheric density model, a 
three-dimensional, static, localized perturbation of electron density characterized 
by a f ew-thous and-kilometer extension in the horizontal plane and a wave-like density 
structure in its interior. This localized perturbation will be cut through, in the com- 
puter simulation, by the DM-to-CSM baseline in successive orbits. We mil attempt 
to reconstruct the TID position and structure by inverting the DM-to-CSM differential 
doppler data obtained with the computer simulation. 

Similarly, we will determine by computer simulation the effect on the DM-to-CSM 
link of turbulent ionospheric regions, such as the aurora oval (to which the radio path 
becomes periodically tangent) and the sporadic E belt (which is periodically cut through 
by the DM-to-CSM path). In this way, we expect to be able to derive the position of 
the boundaries of these turbulent regions from the MA-089 data. 
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6. CONCLUSIONS 


Our computer simulation of the ionospheric measurements to be performed by the 
MA-089 doppler-tracking experiment has shown that at all times the expected value of 
the ionospheric-induced doppler shift in the DM-to-CSM path vastly exceeds the doppler 
noise of the instrumentation system (multipath effects included). For the space-to- 
ground paths, this is also generally true, with some rare exceptions occurring at 
times when the radio path crosses regions of minimum ionospheric density. 

We therefore conclude that this analysis confirms the validity of the MA-089 
ionospheric measurements and encourages the formulation of optimistic expectations 
for the outcome of the experiment. 
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2 , 2 2 8 fc - () l 

4.5736+00 


7b 

25.90 

330.00 

b.BbPF+ni 

3.0OEF4O1 

5.742F-0 1 

1 , 179E+01 


R 0 

2 b. go 

3 3 0. 0 0 

8 . 4 3 3 E * 0 2 

2.E81F + 02 

6. 039F + 00 

1 , 240F+ 02 


Rb 

25.90 

3 3 0 . 0 0 

3 , 7 b J.F ♦ 0 3 

' 1 . 1?1F + 03 

3.4R4F+01 

7, ib6E+02 


QO 

25.90 

330.00 

1 . 3 2 OF * 04 

2.755F +03 

1 .2175+0? 

2.499F + 03 


Qb 

2b. 9 0 

33 0. 0 0 

3 . 1 4 F + 0 4 

4, 482 E + 03 

? . R 8 P F: + 0 ? 

5. 9336+03 


1 0 0 

25. go 

3 3 0 ; il 0 

b . 669F * 04 

5 . c 3 6 F + 0 3 

5.205E+0? 

1 .U89E+04 


lob 

25.9 0 

3 3 0.00 

8 » 6** 2F ♦ 04 

8 . 3 9 7 F + 0 3 

7.934E+0? 

1.8 3 OF + 04 


1 1 0 

25.9 0 

330. n 0 

1 .*niF + 0b 

8 . ~ 9 0 F + 0 3 

1 . 1 02F + 03 

2.283F+04 


lib ?4,9t> 3? 0 . 0 11 1 , 3t>VE * Ob 6 , c 64F + 0? 1.25RF+Q3 2,bR4F+04 

1?0 25,90 330,00 1.375F+05 -b.?45p : -04 l.?*?E*03 2.592E+04 

\pb 2b. SO 330,0 0 l..37bE>tf5 -3,4b9F>04 1.262E+03 2 ,bq?F. + 04 

130 25.9.U 330 , on 1.375F+0S -E.M bh-04 1.262F+03 2.592E+04 

13b 2 b , 90 .3.30,00 1.375F+05 -1,600F-04 1.262F+03 2.5926+04 
140 2b. go. 330, 00 .1 , 39) F+Ob l,47?F:+03 1.H0F + 05 1.31PE+0* 
14b 2b, 90 330,00 1 .465F.+ 06 1,480 F + 0 3 1.2926+05 1,2005+04 

l=,0 ?b.90 330. (jo 1 . 9 3 7 F + 0 b 1 , 4 4bpj+ 0 3 1.468F+09 1.UR7F+04 

1 b 2b. go 330.00 1 . 809F + Ob l ,428F; + 03 1 . b39F * OS 9.819E+03 

180 25. go 3 3 0 , 0 U l . 6 8 0 F + 0 5 1 , 4 0 7 F + 0 3 1.802F+05 R , 855E + 03 

18b 2 b, go 3 3 0 , 0 o l.75up*gb 1.384F+03 1.9576+OS 7.9875+03 

170 25. go 330,00 l,MlRF*Ob 1 , 3 b R (• +03 2 . 1 02F + Ob 7 >229E + 03 

17b; 2 b. go 3 3 0,00 T , RR8F + l.'b 1 . " 2 9 F + 0 ?, 2,2375+05 6.5RQF+03 

lftU 2 b. go 3 30, f)0 l . 9 9 i. F ♦ 0 b 1.2985+03 2 3.6 OF + 0.5 6.079E+03 


1 Rb 

25.9 0 

330.00 

2 , ' 1 5 F + 0 5 

1 . 264 F * 0.3 

2,470 F + 0 5 

9. 711F+03 

1 QO 

25,90 

33 n . 0 0 

? . (* 72 F ♦ Ob 

. 1.2286 + 03 

2.566 F + Ob 

b ,4 q5E ♦ 03 

1 9b 

25.9 0 

33 0. 0 0 

2.1 3HF+05 

1 . 1 R6F+03 

2.647E+05 

5.442F+03 

2 00 

25.9 0 

33 0 . 0 0 

2. 196 F+ 05 

l . 1 4 6 F + 03 

2,712 F+ 05 

5.565F+03 

20b 

25,9 1.!' 

33 0 . 0 >J 

2 .253 F+ 05 

1 .1026 + 0,3 

2 , 7bqF + Ob 

5.8756+03 

2 1 U 

25.9 0 

33 0.11') 

2.307F + 1)5 

1 .nb4F+03 

2.7B8E+05 

6.3R4E+03 

2 1 b r 

25 .9 8 

3.3(1,0 0 

2.35RF+U5 

1 .O04F+03 

2.796F *09 

7.i06E + 03 

2 20 

29.9 0 

330 . on 

2.40 ?F + 05 

9.=n9F+0? 

2.782F+05 

R . 085E+03 

’22bj 

2b . 9 0 

330,00 

2.453F+05 

H,q4RF +0? 

2.74SF+05 

q .244E + 03 

c 3 0 i 

25.90 

330.00 

2.496F+0S 

8 t 7b7E+02 

2.6R4E+09 

1 .06QE+04 

2 3b 

25.90 

330.00 

2.537F+05 

7.7356+02 

2.595F + 05 

1 .24 IE + 04 

’2 4 0 

25.9 0 

33 0 , 1)0 

2.57QF+05 

7 . f • R 0 F ♦ 0 2 

2 , 4 7 R fc ♦ 05 

1.4416+04 

'24b 

25 . 9 0 

3 3 0. 0 0 

2.6075+05 

6. ?9 It +02 

2,3316+09 

1 ,6736+04 

•2R0 

25.90 

3 3 0. 0 U 

2.638F+05 

b.f 68F +0? 

2. 151F+ 05 

l .937E+04 
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*1 

N1 

06 DR 

5 N. n n l 

I'rsDWL 

4 U 

36,6 3 

715 .no 

1 .A'065-ni 

2.6M5F-02 

-2. 152F-03 

1 . 7758-02 

4 b 

3 6,63 

315. 0 0 

4. OH3F-01 

8. 1678-02 

-5.2588-03 

4 . 337E-02 

S ( ' 

36,6; 

3 

115.00 

i.?Pi.F*on 

3.3948-01 

-1 .65PF-02 

1 • 3A7F-0 1 

8b 

36.63 

315.00 

6. } 7 2 F + 0 n 

2,3545*00 

-8.11 OF - 02 

A.6P9F— Ol 

a<> 

36,6. 

3 

31^.00 

5 , 2 4 6 F + 0 1 

2 . C 7 1 F ♦ 0 1 

-6.786F-01 

5 . b 9 7 E + 00 

Ab 

36,6 

3 

315.0 J 

1 • 0 5 n F ♦ 0 2 

2 . , f. q m e' - n 2 

-l .3418 400 

1 . -lOAF + 01 

7 0 

36.63 

315.00 

5 , 4 5 2 (r ♦ 0 1 

- 1 , 9 11 F ♦ n ! 

-7. 5018-01 

6. 1R7F+00 

7b 

36.6 3 

315. 0 0 

6 . n 0 1 F + 0 1 

Z.P17F+01 

-1 .879F + 00 

l .5506 + 0 1 

p U 

3 h . 6 

;•} 

315.00 

6 . 04HF’ + 02 

2.501K +02 

-1.9738+01 

l.b?7fc+02 

nb 

36.6 

3 

315.00 

3.528F+03 

l .0538+03 

-1 . 13RF + 02 

9,3876+02 

qo 

36.6 

3 

315. 0 0 

1 ,p4lF+04 

2 . c 9 0 p ♦ 0 3 

-3.9758+02 

3. 279E+ 03 

qb 

36.6 

!3 

3l5,,.i 'J 

2.O54F+04 

4 , 1 9 4 F + 03 

-9.4368*0? 

7. 7836+03 

loo 

3 6.63 

315. 0 n 

5.329F+U4 

5.305F +03 

- 1 . 7 n l e + 0 3 

1 .4036+04 

lob 

36.63 

315. 0 ) 

5, 1 2 4 F ♦ 0 4 

6. v 1 4 F + 03 

-2.592E+03 

2. 1386 + 04 

3 1 u 

36,6. 

0 

315.00 

) . 1 29F+05 

6.0078 +03 

-3.599F+03 

2.9686+04 

lib 

3 6,63 

3!5.ii 0 

1 • 2 6 76 + 05 

6.1 718+02 

-A . 1 11 E + 0 3 

3.3906+04 

1?0 

36,6 

3 

315. n u 

1 ,2 t} 2F*05 

-5 . 7 14K-04 

-4 . 122E + 03 

3.4006+04 

12b 

36.6 i 

3 1 b . n j 

1 . 0q2F*0 t ' 

-3.3758-04 

-4.122F+03 

3.4006+04 

1 3 0 

3 6. 63 

315,03 

1 .OOiF + nF 

1 ,6RMfc +03 

1 • 7<»4F *04 

3.0786+04 

nb 

36.613 

3 1 5 . ‘) u 

1 ,q99F+05 

1 .555F+03 

4.302F+04 

2.M146+04 

140 

36 ,6|3 

1 1 5 . . j J 

1 .49bF+os 

1 .5366+03 

4.766F+04 

2 .HP06 + 04 

1 4 b 

36 .63 

315.') .1 

1 • 5 9 2 F" + u 5 

1.5036+03 

5.2346+04 

2.549E+04 

150 

36.6 3 

315. 0 0 

1.68/6+05 

.1. .F74F + 0 3 

5.695F+04 

3.0236+04 

15b 

36.63 

3 15, flO 

1 , 779F + 05 

I.F30F +03 

6, 140E + 04 

3. 1026 + 04 

1 A 0 

36 ,6j3 

315.03 

1 .6 70F + 05 

1.7 7.38+03 

6 • 558F + 04 

3.1886+04 

1 Ab 

36 . A 

3 

315.00 

1 . 95 7F + 05 

l', 7 178' + 03 

6.936F+04 

3. H»?F+ 04 

170 

36.6 3 

3 15,0 0 

2. 04 IF* 05 

1 .A 468 +03 

7 . ?6 OF + 04 

3.3866+04 

17b 

36 . 6 

3 

3 1 5 . U j 

2.1 2 F + 0 5 

1 , c 66F + 0 3 

7 . 6 1 IE + 0 4 

3.5036+04 

IPO 

36.6 3 

315. 0 3 

2. }9hf+0 c 

1.4758+03 

7 .673E + 04 

3.6336+04 

1 P.b 

3 6 , 6. 3 

315.00 

2,269^+05 

1 , 3 7 4 F ♦ 0 3 

7.722F+04 

3.7806*04 

IqO 

36,6 3 

315.00 

2.3356+08 

l ,"6lF + 03 

7.6338 +04 

3.9466+04 

iqb 

36,63 

315.03 

2 . 395F + 05 

1.1368+03 

7.376F+04 

4 .1346 + 04 

■? no 

36,63 

315. no 

? . 4 4 M F ♦ 0 5 

9,c,69f +02 

6.913F+04 

4.3476+04 

’c nb 

36.6 3 

315.00 

2.494F+05 

3.4258*0? 

6.2008*04 

4.5928+04 

2 1 0 

36,6 

3 

315. 0 0 

2.532F+05 

6.709F+02 

5.1848*04 

4 . 87 38 + 04 

21b 

36,83 

315,00 

2 , 56 1 F * 0 5 

4.7956+0? 

3.7948+04 

5. 1976 + 04 

•2 ?0 

36 . 6 3 

3 1 5.00 

2.5MUF+05 

2, A52F + 02 

1 .9438 + 04 

5.5756+04 

2 ?b 

36,63 

3 1 5 .0 0 

2 .587F+ 05 

7.0398+02 

3.3396+05 

1 .5266 + 03 

23 0 

36.6 3 

315.0 J 

2.M69F+05 

6. 1 1)76 + 03 

4.8738+05 

-4 . 45 ?F + 0 4 

2 3b 

36.63 

315. 0 0 

3.1 6d F + 05 

5.P628 +03 

2.P1 2F+05 

-3. U 04 6 +04 

24 0 

36,6 3 

315. ) 0 

3 . 4 5 4 F + 0 5 

5 , u 5 0 8" + 0 3 

7.7748*04 

- 1 .4806+04 

24b 

36.63 

315.0 ■> 

3. 722F + 05 

5,1 748+03 

-1 .2086 + 05 

1 *0706+03 

250 

36.63 

315.no 

3 . y 7 OF + 05 

4. 7 32 F+ 03 

-3.1278+05 

1.7466+04 
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M. 
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N 

df or ondm 

dnhwl 

40 

44 .8 6 

3 n o . n 0 

1 . b 3 1 F - 0 1 

2, c 96F-02 -2 • 06 8 E -03 

2.197F-02 

4 b 

4 4 , P b 

300 , oo 

3.P48F-01 

7.PP9P-02 -5,0538-03 

5,3675-02 

«=, 0 

44.8b 

3nfi,uO 

1 .?38Ff + 

3 , ? 8 3 1 - 0 l - 1 • 5 9 3 F - 0 ? 

1 .69PE-01 

5b 

44. Pb 

3 0 n , o'i 

s.QbSF+on 

2,?7f F + 00 - 7 . 7 9 3F ~ 0 2 

8 ,2786-0 1 

6U 

44 .fib 

3 0 0, ijj J 

8 . 0 7 4 F ♦ 0 1 

2.002E+01 -6.52 1 £ • 0 l 

f ,9?PF+00 

6b 

44. fib 

3 0 0. it! J 

1 . U5F + 0? 

2.4826-02 - 1 , 2 8 9 E + 0 D 

1.36PE+01 

70 

44. fib 

300.^0 

S.263F+01 

-1.F568+01 -7.17SE-01 

7 , b ? 6 E + 0 0 

7b 

4 4.8b 

30 0 . 0 -I 

S,B4HF+*)1 

2.5808 + 0 T -1.732E + 00 

1 .840E+01 

R 0 

44. fib 

3 0 0 , -Jj ) 

5 . 5 6 3 F + 0 2 

2.3 008 + 0? -1 .B15F + 01 

1 ,9285+02 

Bb 

44 . Pb 

3 0 0 . 0 0 

3,?4bF+D? 

9,b«8F+0? -1.047F+02 

1 . 1 1 ?E + 03 

90 

4 4.8b 

300 , 00 

1 , 1 4 1 F + (1 4 

2,3821- +03 -3 , b 5 7 F ♦ 0 ? 

3.8R5E+03 

Qb 

44 ,Pb 

300,00 

2 . 7 l 7 F ♦ l) 4 

3 , 8 5 7 F + 0 3 -R.bP2F + 0P 

9.223E+03 

1 no 

44,8b 

300.00 

4 . 9 0 1 F + 0 4 

4.7RfF>n3 - ] , 5b5E +03 

1,6625+ 04 

1 ob 

44 . fib 

3 0 0 . o o 

7, 4 7 2F + 04 

b , c 3 0 h* + 0 3 -2.385E + 03 

?.b348. + 04 

1 1 u 

44 , P 6 

300,00 

1 . 0 39 F ♦ OB 

5 . t. ? 4 F + 0 3 - 3 • 3 1 1 E + 0 3 

3.bl7F+04 

1Tb 

44 . Pb 

3 n o . u o 

] . ] 8 JF ♦OF 

5 . f 7 5 F + 0 2 -3.782E + 03 

4 . U 1 RE+ 04 

1 9.0 

44 .fib 

300,00 

1 , 189F+05 

-5.043F-04 - 3 . 7 9 3 E + 03 

4.U2QE+04 

1 ?b 

44 .Pb 

3 0 0. 0 0 

1 . 1 89F + Ob 

-3.?b4f-04 -3.793E+03 

4 , U29F + 04 

170 

4 4.86 

300.00 

1 . ?. 5 4 F ♦ 0 5 

1 . 9 P 5 F * 0 3 1.4P3F + 04 

4, I25E+04 

17b 

44.8b 

300 , JO 

l . 3S3F+05 

1.9796+03 1.531E+04 

4.2POE+04 

140 

44 . Hb 

3 0 0. 0 ,3 

1 .45 lF + 05 

1 , 9 6 3 F + 0 3 1 .5R8F + 04 

4 , 4 4 ? E ♦ 0 4 

14 b 

44.8b 

300.00 

1 .349F+0 C 

1 . 9 3 5 F + 03 l.bblF + 04 

4.61 1E+04 

mo 1 

44 .fib 

3 o n , u u 

1 . 645F+05 

l , 8 97 F +03 1.7486+04 

4 , f 8 b F. + 0 4 

1 5b 

44,8b 

30 0.00 

1 .738F+05 

1.P47F + 03 1 ,H44fc + 04 

4.970E+04 

If U 

4 4,8b 

3 0 0. 1) 0 

1 § 8 ? v F * 0 5 

1.7B4F+03 1.9448*04 

5. 160E+04 

lfb 

44.6b 

300,00 

.1 ,P16F*0 C 

1 *7 03F + 03 2.040F + O4 

5.357E+04 

170 

44 .Pb 

36 0. ij 0 

2 . (i OOF * 08 

1 ,b 1 88 + 03 2. 123E + 04 

5 . bf ] E ♦ 04 

17b 

44,86 

3 0 0 . O'* 

2 . 0 7 bF + 05 

1,5128 + 03 2*181 F + 04 

5.7705*04 

1B0 

44,8b 

300,0 J 

.} 5 lF + 06 

l .7 8PF + 03 2 • 1 9 8 F + 0 4 

5 . 9 8 3 E + 0 4 

1 flb 

44.6b 

300. 00 

2.217F+05 

1 # 2 4 8 F + Q 3 2.1538 + 04 

6. iPRE + 04 

1 <5 0 

44.86 

300,00 

2,?7bF+05 

l .0798+03 2 . 0 1 9 F * 0 4 

6 .4 1 4E' + 04 

1 Pb 

44.8b 

3 0 0.0 ) 

? . 3 24F + OS 

0.P55F+02 l.tbOF+04 

6.6 26E + 04 

c n o 

44 .Pb 

30 0 . 0 0 

? , 36 3 F + 0 5 

b . F P 6 f + 0 ? 1.3248 + 04 

6.8305*04 

c ob 

4 4 . 6 b 

3 0 0.0 

2. 3B9F+05 

b. 93 IF + 02 6. 3888 + 03 

7.0?OE+04 

c 1 0 

44.8b 

3 0 0 . n 0 

P . 4 () 1 F ♦ 0 3 

7.424F+01 -4.063F+03 

7.1856+04 

2 1b 

44 ,P6 

300,00 

2.522F+05 

3.0018+03 5.912F+04 

5,0688+04 

2?0 

44 .Pb 

3 0 0 , ,t U 

? ,bb9F + 05 

2.PP8F+03 -3.41 4E +03 

3 .6P6E + 04 

■2?b 

44. Pb 

300.00 

?. 8 UP + 08 

2.7678+03 -6.4218+04 

2. 7P7F+04 

c 7 0 

4 4.8b 

300,00 

P.946F+05 

2,6 048 + 03 -1.227E + 05 

1 , 7P7E + 04 

2 7b 

44.8b 

3 n () . it 0 

3 . { 1 7 1 F + O' 5 

2 .4 05p + 03 -1 .7P5E + 05 

fi ,H01 E + 03 

2 4 0 

44.8b 

30 0 .00 

3 , 1 8bc + O t ~ 

2, 1 6 3 F’ ♦ 0 3 -2.31 2F + 05 

8 , 9 44E + 02 

2 4 b 

44 .Pb 

3 0 0 . 0 0 

3 . fi 8 b F + 0 5 

1 .8 70 F +03 8068 + 05 

-5,7585*03 

'250 

44,6b 

300.00 

3 , 3 7 1 F + 0 5 

1,514 F~ + 03 -3.2708*05 

-1.111E+04 
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A 0 

bO . 0 J 

POft.il i 

1 ,b5 7F-:t)l 

2. *,946-02 1.631E-03 

P.O91E-02 

4b 

b i) , n j 

?«{=>.:»•) 

3. 793F-0T 

7.C07E -np 3 . 9 HAF -03 

6.3736-02 

8 0 

ftn.yj 

P Oh . H 

1 . ! 098*00 

3,1536-01 1.PS6E-0P 

2. U73E-0 1 

8b 

bn.n.'t 

POb.tr! 

b . 733F+0 0 

2.1876+00 6.1A5F-0? 

1 .0146+00 

6 0 

b 0 . 0 J 

P Pb „ "1 1 

a , o 7,lF + n } 

1 , 9 ? 3 6 + 0 ] ft.lAPF-O] 

8.4P4F+00 

6ft 

0 , 0 3 

ppb.no 

9. 74 9 F + 41 

2.2416-02 1.016F+00 

1.6766*01 

7U 

b 0 , 0 3 

? « 5 . -j 0 

ft . OAbF+01 

-1.7916+01 ft. 6 366- 0] 

9,2996 *0-0 

7 b 

b n , n 3 

Pob.no 

b . OAbF* 0 1 

2. 3 17 F + 01 1.3 OAp* On 

2 . I ft 1 F * 0 1 

ftO 

SO. D8 

Poft . o J 

ft, 0 2bF* OP 

P. 0776 +02 1.3646*01 

2 , 2ft0F+ 02 

P b 

bO .08 

?o b . o ;i 

P.9 3 OF ♦on 

0 . 7 A 9 (- +02 7.0666 + 0 1 

1 . 298F* 03 

90 

b 0 , 0 3 

P 8 ft , ,)0 

1. , o 3 0 6 + 0 A 

P. 1ft It +03 2. 7406 + OP 

A .3336*03 

9 b 

b 0 . 0 3 

POb. on 

P , A ft A F ♦ 0 A 

3.4H4F+03 6.5226 + 0 ? 

1.0766+04 

1 nO 

-» 0 . 0 3 

?ob . no 

A , A 2 6 F + 0 A 

4,3236+03 1.175E+03 

1 .9396*04 

l nb 

b 0 , 0 3 

POb . i 

6 • 7 A H F + 0 A 

4.9956+03 1.7926 + 03 

2.9866*04 

J 1 u 

b (1 . 0 3 

p p b . ( i n 

9 , ^806 + 1)4 

4.9R9F+03 2.4076*03 

A • 1 046*04 

Ho 

bn . C 3 

poo , «:} o 

l . 0 6 9 F ♦ 0 ft 

b . 1 ? 5 F ♦ 0 P ?. 0416 + 03 

4.9886+04 

1 ?0 

bO.0 3 

? 0 b • <J 0 

) . 0 7 J F + 05 

-4.0446-04 2 • 8 4 9 F + 03 

4. 70J6+04 

1 ?b 

bO , 03 

pob.a ) 

l . 0 7 j F + 0 ft 

-3,1 3 b F - 0 4 ?.849F*03 

4 • 7 0 1 F + 04 

no 

b 0 « 0 J 

? w S . 0 0 

1 , ] 8PF+0R 

1 .0476+03 -5.961P+0? 

ft. 3786+04 

nb 

b 0 j, fi 3 

POb.M J 

l .PP46+05 

1 .P39F +0,3 - 3 * b 7 7 £ + 0 3 

8.0426+04 

l 14U 

b 0 . 0 3 

POb. 1 J 

1 , 3 1 6 F ♦ 0 ft 

1 . P 226 +03 -6.3676 *03 

ft.9lftF^04 

, lAb 

.1 ft o 

b 0 » 0 3 

p oft. no 

l .406F+05 

1.7946 +03 -R. 9196 + 03 

6. 1 9 7 F ♦ 0 4 

ft 0 . 0 3 

?»ft , n 0 

1 ,49bF+0b 

1.7866+03 -1.119F+04 

6,4886*04 

nb 

b () .03 

POb . 0 0 

1 , 5 H 2 6‘ + 0 ft 

1,7076 +03 -1.3126 + 04 

6.7786+04 

no 

b 0 . C 3 

POb. (j J 

1 , 66b F + Ob 

1,6476+03 -1.4686+04 

7.0726*04 

16b 

bn , 0 3 

POb . u it 

1 , 7 A h F * 0 ft 
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